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In this methodological note we consider, in a weak-field limit, a relativistic linear motion of 
two particles with opposite signs of masses and a small difference between their absolute values 
mi, 2 = ±(/i ± A/i), n > 0, |A^| <C fi. In 1957 H. Bondi showed both in framework of Newtonian 
analysis and in General Relativity that when the relative motion of particles is absent such a pair 
can be accelerated indefinitely. We generalise results of his paper to account for a small nonzero 
difference between velocities of the particles. 

Assuming that the weak-field limit holds and the dynamical system is conservative an elementary 
treatment of the problem based on the laws of energy and momentum conservation shows that the 
system can be accelerated indefinitely, or attain very large asymptotic values of the Lorentz factor 
04 ' 7. The system experiences indefinite acceleration when its energy-momentum vector is null and the 

mass difference < 0. When modulus of the square of the norm of the energy-momentum vector, 

■ l-^'^li is sufficiently small the system can be accelerated to very large 7 oc |A''^|~^. 
It is stressed that when only leading terms in the ratio of a characteristic gravitational radius to 

f~| , the distance between the particles are retained our elementary analysis leads to equations of motion 

■ equivalent to those derived from relativistic weak-field equations of motion of Havas and Goldberg 
1962. 

I Thus, in the weak-field approximation, it is possible to bring the system to the state with extremely 

high values of 7. The positive energy carried by the particle with positive mass may be conveyed to 
other physical bodies say, by intercepting this particle with a target. Suppose that there is a process 
^ , ' of production of such pairs and the particles with positive mass are intercepted while the negative 

meiss particles are expelled from the region of space occupied by physical bodies of interest. This 
$JL( ' scheme could provide a persistent transfer of positive energy to the bodies, which may be classified 

as a 'Perpetuum Motion of Third Kind'. 

Additionally, we critically evaluate some recent claims on the problem. 
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I. INTRODUCTION 



Bondi 1957 [1] pointed out that in the Newtonian approximation two particles with opposite signs of masses at 
, rest with respect to each other accelerate indefinitely in an inertial frame. This process is allowed by the laws of 
OJ ■ conservation since the kinetic energy and angular momentum of such a system are conserved being exactly zero while 
the potential energy depends only on relative distance between the particles. In the same paper he generalised this 
result by finding an appropriate static accelerated solution in General Relativity and discovered that a uniformly 
accelerated pair of particles with opposite signs of masses must have a mass difference determined by the fact that 
constant in time particle accelerations must be different to keep them static with respect to each other. 

It is trivial to show that in the Newtonian approximation (see the next Section ) when the two particles with 
opposite signs of masses have a relative velocity, its value is approximately conserved. As a result, the acceleration 
period is finite and the pair as a whole being initially at rest gains a finite value of velocity. We also show that when 
initial relative velocity of the particles is sufficiently small the pair can be accelerated to a relativistic speed. 

In the next Section we consider the problem in the relativistic setting and generalise Bondi's analysis considering 
pairs of particles with opposite masses and small difference between their absolute values: mi 2 = ±(/i± A/i), > 0, 
|A/i| <C ^ having an initial relative velocity Vm in a fixed lab frame where the pair as a whole is initially at rest. 
We assume that gravitational interaction is weak and, therefore, G/i/(c^i5i„) <C 1, where Din is an initial separation 
distance between the particles. Also, for simplicity, in the relativistic treatment, it is assumed that the orbital angular 
momentum of the system is equal to zero and the motion is linear. 

We analyse this situation by elementary means. The equations of motion are obtained from the laws of energy and 
momentum conservation. It is assumed that energy and momentum of the system in a Lorentz frame instantaneously 
comoving with the motion of the pair are given by the Newtonian expressions and that they form time and spacial 
components of a local four-vector. Then this energy-momentum vector is projected onto the lab frame. Since energy 
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and momentum in the lab frame are conserved under the assumption that gravitational radiation from the system is 
insignificant we get two equations of first order in time fully describing dynamics of the system. We also show how 
to derive an equivalent pair of second order equations considering the Newton's law of gravity in a frame accelerating 
with the particles. 

It is shown that the pair as a whole always have a positive acceleration with its asymptotic value being either zero 
or a nonzero constant depending on initial conditions. The relative distance between the particles can either have a 
turning point or increase monotonically. The system accelerates indefinitely when the mass difference A/x < and the 

norm of the energy- momentum vector N — (2A^c^ + ^jy—Y ~ l^^'^in — Sijad, accordingly, the energy-momentum 

vector is null. In this case the relative distance increases monotonically. When is sufficiently small, for the initial 
conditions corresponding to the monotonic behaviour of the relative distance the acceleration period is finite but the 
asymptotic value of the Lorentz gamma factor is large being proportional to |A'^^|~^. 

Such pairs can play a role in realisation of a hypothetical effect, which we called 'Perpetuum Motion of Third Kind' 
[2], hereafter PMT. In its most general formulation this effect is a possibility of a persistent energy transfer from a 
subsystem having negative energy to a subsystem with positive energy, in classical theories where negative energy 
subsystems are possible. Indeed, the positive mass particle can, in principle, be used to transfer positive energy to 
other physical bodies after the pair has been accelerated to high values of the Lorentz factor. Iterating this process 
as many times as we need we can extract as much positive energy as we wish. Note, however, that it is not the 
only 'working model' of PMT, and that, in principal, in order to make PMT we need neither systems with negative 
rest mass nor gravitational interactions. As is shown in [2] it suffices to have a medium violating the weak energy 
condition with certain additional properties and mere hydrodynamical interaction 'to construct a PMT'. 

Additionally, we comment on several statements of paper [3] , where the Kepler problem for a binary with opposite 
signs of masses has been considered which may, in our opinion, lead to misunderstanding of the problem. 



II. NEWTONIAN TREATMENT OF THE PROBLEM 



At first let us consider the problem in the Newtonian approximation where mutual gravitational accelerations acting 
on the particles of masses mi and rn2 are given by the conventional expressions: 



Gm2 .. Gmi 



|D|3 



|D|3 



(1) 



where are position vectors of particles with indices i = 1,2 and D = ri — r2. Setting fi = Gmi = —Gm2 we obtain 
from (d]) 



IDP 



D, v = 0, 



(2) 



where V = ^{ii + Y2) and v = D. From equation ^ it follows that when v = at some moment of time it remains 
zero in the course of evolution of the system. Thus, in this case the relative separation distance D does not change 
during the evolution and the system ever accelerates as a whole, with the acceleration vector 



V = 



DP 



D 



(3) 



being constant. The laws of conservation are nonetheless respected since the kinetic energy and momentum of the 
system are precisely zero while the potential energy depends only on the relative separation distance^. 

When vit = 0) = Vi„ 7^ the absolute value of the relative distance changes with time. Accordingly, the absolute 
value of the acceleration changes as well and eventually decays provided that (Di„ • Vi„) ^ — |vi„||Di„p, where 
Di„ = T>{t = Q). We have 



D = Vi„i + 



(4) 



^ Note that it is easy to show that the same motion can be reahsed in a system containing particles provided that the total mass of 
the system M = ^22ii—i n ™' ~ ^ ^^'^ positions of the particles are chosen in a special way. Say, for a system containing three particles 
their relative positions must form an equilateral triangle. 

^ Clearly, when (Din • Vi„) = — |vi„llDi„| the particles collide. 
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and, thus, integrating equation ([3]) we obtain 



1 u fi/^ + ar 
|V(t)| = , J2(l ^), (5) 

where we use the dimensionless time r = J-Q^t, Din = |Dm|, '^m = |vi„|, e = j./'^a , a = (v™ • Di„)/(i;j„ An), and 



A = Ve + t2 + 2aei/2r. 

Note that when the system moves along a straight hne with increasing vahie of |D|, and, accordingly, a — \ equation 
© yields 

In the limit r — )■ oo we get from equations ([5]) and ([6l) 



^ |V(r ^ oo)| = ^/-_— (7) 

" 1 + a DinVin 



It follows from ([7]) that when 



«»<.c„* = \/^;£^, (8) 



the asymptotic value of velocity of the system, Vaoi formally exceeds the speed of light, c. Clearly, a relativistic 
approach to the problem is to be used in this situation. 



III. RELATIVISTIC TREATMENT 
A. Derivation of dynamical equations 



In order to keep our study as simple as possible let us consider in the relativistic case only the motion along a 
straight line with increasing value of D (a = 1). Additionally, in this Section we use the natural units setting the 
speed of light and the gravitational constant to unity. However, unlike the Newtonian case, here we would like to 
consider particles having a small mass difference: mi. 2 = /i± A/x, where it is assumed below that /i > and |A/i| <C /i. 

It is useful to introduce two local frames and the respective coordinate systems 1) a fixed lab frame with global 
Lorentzian coordinates (x, t) and 2) a local Lorentzian frame instantaneously comoving with the motion of the point 
R{t) — ^{xi{t) + X2(i)), where xi{t) and X2{t) are positions of the particles in the lab frame, with associated 
Lorentzian coordinates (a;^°™, t'^"™). Is is assumed that at some particular moment of time t — t^, coordinates of the 
event (t*, R{t = <:»)) in the comoving coordinate system are equal to (r, 0), where t is the proper time associated with 
the world line (t, R{t)). Hereafter, the world line (t, R{t)) is referred to as "the reference world line". 

When = r the positions of particles are given by a;f°™(i'^°'"), their velocities are w^"™ = dF^^TT- Let us 
also introduce the relative position and velocity in the comoving coordinate system D — a;™™ — a;™™, = D. 
Without loss of generality we assume hereafter > 0. When the relative separation remains sufficiently small 

along the reference world line we have approximately a;™™ = —xf"^. 

In the global coordinates at the time slice t = t^, the velocity of motion of the system as a whole is given by 
^ ~ ^(^•^i ^^2)(i = t*) while the relative position and velocity of the relative motion are Diab = xi{t^) — X2{t^) 
and V = f^D. 

Introducing the Lorentz gamma factor 7 = ^-^^ y'i associated with the reference world line we may write in the 
limit of small separations 

D{t-n = lD,a,. ^=7, (9) 

and, accordingly, 

V^°'^=Jj^hDiab)^J^V+J^Diab. (10) 
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Supposing below that, on one hand, the relative distance D ^ fj,, and, therefore, a weak-field approximation holds 
and, on the other hand, it is not too large for the local Lorentzian coordinates to be adequate and, respectively, 
equations (|9l fT0|) to be valid we can use the Newtonian expression for the energy, Ec, and momentum, Pc, of the 
system in the comoving frame at the time t'^°^^ = r 

E,^2Afi + ^, Pc^fiD, (11) 

where dot stands for differentiation w.r.t. the proper time r. 

In the same limit and Pc represent time and spacial components of a local four vector, and, therefore, they 
values in the lab frame, E and P, respectively, can be obtained from (1111) by the standard Lorentz transformation. 
We have 

2 2 

E = -f{2A^^+^ + V^,D), P = j(^,D + V{2A^^ + !^)), (12) 

where it is assumed that the velocity of the systems as a whole, V, is a function of the proper time r. Since energy 
and momentum in the lab frame are obviously conserved, equations (I12p fully describe the dynamics of our system. 
They should be solved subject to the condition that the system is initially at rest with respect to the lab frame: when 
T = we have V = and 

E = E„, = 2A^ + P^P^n^ (13) 

where Din and Vm are initial separation distance and relative velocity, respectively. It is assumed below that Wi„ > 0. 

Although our derivation of dynamical equations (|12|) may look somewhat heuristic it is worth mentioning that 
when terms next to the leading order in /i are discarded they can be derived from the precise weak-field equations of 
reference [4] in the limit of small separations and |A/i| <C fJ. ■ 

It is convenient to transform equations (|12p to another form using their linear combination E — VP and calculating 
square of the norm of the energy- momentum vector, N'^ = E^ — P^. We get 

E~VP=^-f-\2A^i + ^) (14) 



and 



iV^ = (2AM+^)2-A*^(I))^ (15) 

We also obviously have N'^ = (2A/i+ )^ — M^^^n- Note that contrary to the usual situation the energy-momentum 
vector can be null, time-like or space-like, depending on initial conditions. 

Equations (|12p are first order integrals of two second order in time dynamical equations. One of these equations 
can be obtained from (|15p by differentiating this equation over r with the result 

^ = -^-;d3' (16) 

and the second one by differentiating either of equations and using 

^'y - ^ (17) 

Equations ((TB)) and (|17l) can be obtained from other independent qualitative arguments. The derivation of the 
second order dynamical equations, which relate dynamical variables with different values of time coordinates is not, 
however, convenient in the local Lorentzian coordinates introduced above since these coordinates are defined with 
respect to some particular event on the reference world line and, therefore, the definition is different for different 
events along this world line. It is much more convenient to use a coordinate system, where the proper time r plays the 
role of coordinate time. For that let us consider another, the so-called local Fermi- Walker coordinate system (r, y), 
see e.g. [5], where the proper time r is the coordinate time and the unit vector in the spacial direction y is always 
perpendicular to the four velocity along the reference world line. The coordinates of the reference world line in this 
coordinate system are simply (r, 0). 

At the time slice f^"™ = r the local Lorentz coordinates and the Fermi- Walker coordinates coincide: xf™ = yi,2, 
but the Fermi- Walker coordinate system is accelerating with respect to the local Lorentz coordinate system with an 
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acceleration g{y). Clearly, g{y = 0) must coincide with modulus of four acceleration of the reference world line with 
respect to the lab frame. The equations of motion in the Fermi- Walker coordinates are assumed to be determined by 
the Newton's law ([Ij with added acceleration term —g, which accounts for the fact that this system is not inertial: 

T A/i 

yi,2 = -p^2 9, (18) 

where Dpw = yi — 1/2 and we take into account that the acceleration term depends on the coordinate y, e.g. [5]: 
g = a + a^y. For the average distance Y = (j/i +j/2)/2 to be at rest Y{t) = the acceleration term a must be balanced 
by the gravity term „^ . We get 

FW 

(19) 



^FW 



Taking into account that in the lab frame the spacial coordinate of four acceleration is related to a as = 7a we 
have 

U^ = j'V = ja, V=-4^- (20) 

7 ^FW 

It is clear that this coincides with equation pT|) . 

The dynamical equation for the relative distance Dpw directly follows from and <\1\)\ : 



Vpw - -—^ a Ufw-- j:;2 7^3 — • i^l-' 

^FW ^FW ^FW 



It coincides with (jl6p . 

The last term on the right hand side of (j2ip is due to the non-uniform acceleration force appearing in the Fermi- 
Walker coordinates. Because it is c>c /i^, technically, it is a post- Newtonian term. Since we consider the gravitational 
force in the Newtonian approximation in (|21[) it is important to check whether or not post-Newtonian corrections 
to the gravitational force are comparable with the acceleration term in (j2ip . In fact, as is described in standard 
handbooks, e.g. [6], the post-Newtonian corrections are either proportional to A/i or j/i 2. The mass difference and 
velocities are assumed to be small and therefore, the terms in (|2ip arising from the post-Newtonian corrections appear 
to be small compared to the terms taken into account. 

From pTjl it follows that when the mass difference is negative and Dpw = 2|A/i| the particles are at rest with 
respect to each other. In this case the Fermi- Walker coordinate system locally coincide with the Rindler one and 
the particles accelerate indefinitely. Thus, unlike the Newtonian case considered in the previous Section the particles 
accelerating indefinitely being at rest with respect to each other must have the small mass difference. This effect was 
first noted by Bondi 1957 [1]. It is obviously due to the non- uniform character of the acceleration term. 



B. Solution of dynamical equations 

Since equation contains only V and D it can be used to express V in terms of D 

EP , Er 



y=EfTl^(^T^.M^^), (22) 

where Ec is expressed through D in equation (|11|) . E and P are given in equation (|13|) . As we discussed above we 
assume that at the initial moment of time t — t = we have V — 0. That means that initially we have to choose the 
sign (— ) in (|22|) . However, under certain conditions discussed below the direction of motion of the particles relative 
to each other and, accordingly, D, changes sign. At the turning point D = we have = E'^. Since velocity V 
must grow monotonically according to pT| we must take the sign (+) in (1^ after the turning point. 

On the other hand equation (|15p contains only D and its derivative with respect to the time r, and, therefore, it 
can be integrated to obtain the dependence of D on time. Explicitly we have 

^ xdx 



fuf-^ ^/m, (23) 



where 



R{x) = (^2 2A^2;)2 - N'^x^. 



(24) 
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FIG. 1: The dependence of the potential U on the spacial coordinate D. The sohd curve corresponds to the case A/i > while 
the dashed one to the case A/i < 0. 



The integral in p3|) can be evaluated by a standard substitution to give an explicit relation between r and D. 
However, the final expressions are rather cumbersome and we do not show them here. Instead, in general, we analyse 
qualitatively solutions to (|15|) based on analogy between this equation and an equation describing a motion of a 
particle in a potential well. For that we bring (|15|) to a standard form 



where 



^ + U{D)=£^ U{D) = -^-^, (25) 

2^2 "2 A„ 2Dl- ^'^^ 

Introducing natural units C/ = 'KJIP^ ^^'^ ^ ^ ^^^-^ '^^ express U in terms of in a very simple form: 
U = =F-^ — where the sign — (+) corresponds to > (A^ < 0). The dependence U{D) is shown in Fig. 1. 

At first let us consider in detail an important case of zero norm of the energy-momentum vector, = 0, and set, 
accordingly, P = E. A simple analysis of equation (|22p shows that in this case there are no turning points, the relative 
separation D grows with time and the value of y = 1 can be achieved in the asymptotic limit r — )■ cxd. Therefore, in 
this case the system may accelerate indefinitely. 

When = equation ([^ simplifies to 

E^-E,\E,] {E^ + E,\E,\) 
^= E^+E-. ' 2EE. (^^^ 



and from equation we get 



^ ,{2Af,{D - Dra^n) - t^^logi^-^-^^^)). (28) 
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From equation ((27|) it follows that when Ec > the indefinite acceleration is possible only if Ec ^ when r — )• cjo 
and from the expression for E^. pT|) it is seen that the mass difference A/i must be negative for that. We consider 
below only this case in detail. When |A/i| ^ Ec ^ provided that D Dcrit — M^/(2|A/i|). Equation ([28]) tells 
that the logarithm on the right hand side diverges when D — ^ Dent- That means that this limit does correspond to 
the limit t — )> cx). Let us estimate the dependence of the Lorentz factor 7 on time in this case. To do so, we introduce 
a new variable A = Dent ~ D and substitute it to (I28|) assuming that it is small. We get 



4A^2 2|Am|A 
and, substituting this result into equation (j27p we have 



log( "';,y,y'" ). (29) 



4|A/i|A„ '''''' ~^ 



7 « ,1^^,^ exp — y-T. (30) 



Equation (j30p tells that when D Dcrit acceleration is exponentially fast. 

The degenerate case A/i = must be analysed separately. In this case from ([^5]) we have 



2n 

and the distance D increases indefinitely with time. From equation (I27[) we obtain 



l^{D^^Dl,J, (31) 



7«^7^V27. (32) 
Now let US turn to the general case A^^ 0. Setting _D = in f[25|) we get a general equation for the turning points 




Equation p3p tells that the turning points exist only when N'^ > 0. Their number depends on signs of £ and A/i. 
When A/i > the potential U{D) is negative, see Fig. 1, and therefore, the relative motion is finite for £ < with 
one turning point ^ 

In the opposite case £ > 0, and, accordingly, N < 2A/i, the motion is unbound and the relative distance D grows 
indefinitely with time. 

2 

When A/i < the potential U{D) acquires positive values for D > 4]^^, see Fig. 1. It tends to zero when D ^ 00 
and has a maximum at Z? = Dcrit- Note that from the condition U{Dcrit) = £ = we get iV^ = there. The 

character of the relative motion depends on whether £ is negative, belongs to the interval < £ < ^^r- corresponding 



to < < 2|A/i|, or £ > — ^ and, accordingly, N < 0. When the energy £ is negative the motion is bound with 
one turning pomt 

In the intermediate region < £ < "^^1^ there are two turning points 



Let us remember that we consider only positive values of D. 
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When Din < -D_ the motion is bound while for Din > D grows indefinitely. Finally, when iV^ < the motion is 
always unbound. 

When the motion is bound the velocity D changes sign after the turning point and in this case we should use the 
sign (+) in (1^^ . Taking into account that D is decreasing after the turning point and that Ec oc D~^ we see from ([^^ 
that the velocity V 1. The particles tend to collide. However, our assumption that D ^ fi breaks down in this case 
and we cannot describe the motion at scales D ^ ^ within the framework of our formalism. Note that we consider 
in this study only pairs of particles with strictly zero angular momentum. In the situation when the particles have 
a small but nonzero angular moment they would miss each other and after certain moment of time the distance D 
would become negative. In this case the analysis of this paper can be repeated without any major change for negative 
values of D and one would conclude that for such parameters of motion there is another symmetric turning point at 
negative values of D. Thus, the relative motion of a pair of particles with small but nonzero orbital momentum would 
be periodic much similar to the case of ordinary particles with positive masses. 

Now let us consider the case of the unbound motion and estimate the maximal value of the Lorentz gamma factor 
the system can reach. As follows from our previous discussion when N'^ ^ the distance D grows indefinitely. That 
means that the energy in the comoving frame, Ec, must tend asymptotically to 2A/i. Note that when A/i < the 
asymptotic value of Ec is negative. We have from setting Ec = 2Afi there 

^ ^ 4A^2 ^ ^2 (gy^^^- 2AmV4A,.^ - N^). (37) 

Equation (1571) tells that when Afj, > the last term in the brackets is negative and the asymptotic value of velocity 
is smaller than 1. Large values of V can be achieved in the opposite case Afi < assuming |A^^| ^ AfP. In this case 
we expand expressions in ([57]) in the Taylor series in \N'^\/ Afi^ to obtain 



and, accordingly 



, 1 -4|AMi^l (39) 



Equation p9p tells that for fixed values of E and A/it < the gamma factor can be made arbitrary large by choosing 
arbitrary small values of \N'^\. This conclusion is in agreement with our previous finding that the system accelerates 
indefinitely when N'^ = 0. 



IV. METHODOLOGICAL COMMENTS 



Here I would like to make comments on several methodological issues related to the problem. 

1) At first glance the fact that the 'average' position of the pair (xi + X2)/2 always grows with time may seem 
to be in contradiction with the law of conservation of the centre of mass of the system. This contradiction is 
resolved by observation that for the system containing particles of opposite masses position of the centre of mass, 
_R, is determined by a difference of positions of particular particles. Say, in the Newtonian approximation we have 
R — mixi + 1112X2 = {fj, + A^)xi — (/i — AiJ,)x2- In the relativistic case the situation is analogous for systems with 
N'^ > 0. In the opposite case the notion of centre of masses is ill defined. Indeed, introducing the velocity of a 
coordinate system, where the centre of mass is at rest in a standard way as Vcm = P/E [6] we see that when iV^ = 
Vcm = 1 and when iV^ < Vcm formally exceeds the speed of light. It is obvious that the notion of the centre of mass 
is redundant in both cases. 

2) In Introduction of their paper the authors of [3] claim that the conception of PMT put forward by the author of 
this note is related to the problem of indefinite acceleration of two gravitationally interacting particles. This statement 
needs, in my opinion, a clarification. First, let me note that as it is discussed above even in the case when only a finite 
acceleration of the particles is attained PMT is still possible in a situation where production of such pairs is provided 
by some physical mechanism. Second, the conception of PMT, in general, does not rely on gravitation interactions at 
all. In particular, in paper [2] I consider a model where there is a continuous flow of positive energy from some spacial 
regions having negative energy to other regions with positive energy provided by hydrodynamical effects. In this model 
the space-time is assumed to be flat and gravitational interactions are absent. Moreover, in order to construct a PMT 
it is not necessary to invoke objects having negative rest masses. It is enough to consider a medium with positive 
comoving energy density violating the weak energy condition [2]. Additionally, there are ways of constructing PMT, 
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where gravitational interaction plays a totally different role, say, transferring the energy from a non-stationary system 
having negative mass to gravitational waves, as for example in the model of a rotating relativistic string connected by 
two negative mass monopoles [2], [7]. The effects related to the dynamics of free negative mass particles are clearly 
irrelevant to such systems. 

3) The authors of [3] claim that it is impossible to obtain, in principal, an indefinite acceleration of the system 
containing two particles with opposite signs of masses. One may think that this clearly contradicts to the Bondi's 
result [1] and the results of this paper. The conundrum is resolved by observation that the authors of [3] consider 
only relative motions while Bondi's analysis as well the analysis in this note also deal with the motion of the pair of 
particles as a whole with respect to an inertial frame. 

V. CONCLUSIONS 

In this note we show by elementary means that in the weak limit approximation a pair of particles having opposite 

values of masses can be accelerated indefinitely provided that the energy-momentum vector characterising the system 
is mill. The system can also be accelerated to arbitrary large Lorentz factors when the mass difference A/u < and 
the norm of the energy-momentum vector is sufficiently small. 

Assuming that there is a process of production of such pairs and that the positive mass particles are intercepted 
with a target while the negative mass particles fly away it is possible to transfer to the target any desired amount 
of energy. In a more natural situation one can also consider a theory whore the positive and negative mass particles 
interact differently with a normal matter. A general situation of this kind where there is a persistent transfer of 
energy from a subsystem having negative or almost zero energy (like this pair ) to a subsystem with positive energy 
was dubbed by us 'Perpetuum Motion of Third Kind' (PMT) [2]. Note, however, that it is just a classical analog of 
the well known instability of a quantum system with a number of negative energy states unbound from below. 

The question of whether the existence of PMT or ever accelerating pairs of particles is a paradox depends, in our 
opinion, on definition of what paradox is. On one hand, for example, Bonnor 1989 states 'I regard the runaway (or 
self-accelerating ) motion ... so preposterous that I prefer to rule it out by supposing that inertial mass is all positive 
or all negative' [8]. Clearly, existence of PMT can also be classified as a kind of runaway. On the other hand, no 
laws of physics are broken in such systems. We believe that the existence of runaways of these kinds is dangerous for 
theories where they present. To exemplify, an indefinite concentration of energy of different signs in spatially separated 
regions could lead to a highly inhomogeneous space-time hardly compatible with presence of any life. Therefore, such 
theories should be ruled out though some additional study of them in General Relativity may be of some interest. 

Since in our approximation only linear metric perturbations and one next-to-the-leading order term determined by 
the acceleration of the pair as a whole are taken into account, it is interesting to estimate what kind of corrections can 
be obtained by considering other higher order terms quadratic in metric perturbations? For a non-relativistic motion 
with y <i; c for this purpose one can use the well known Einstein-Infeld-Hoffmann equations of motion (e.g. [6]). In 
this way it is convenient to consider particles with a large mass difference as well as systems with non-zero angular 
momentum. There are, however, many corrections, which are absent in such a treatment, notably the emission of 
gravitational waves. Therefore, a self-consistent relativistic treatment of the problem in the next to the weak field 
approximation must be based on the second-order formalism of Havas and Goldberg 1962. Such an approach is left 
for a possible future work. 

Although in this paper we consider only particles with no internal structure our analysis may also be vahd for a 
pair of extended objects with total energies of opposite signs provided that they have a sufficiently large separation 
distance and that their relative velocities are sufficiently small. For example, Deser and Pirani [9] considered the 
behaviour of systems with all possible incrtial/gravitational mass signs and noted that a pair of geons having opposite 
signs of their total energies would behave as a pair of point particles in the corresponding limit. 

It is also interesting to point out that the notion of 'Perpetuum Motion of Third Kind' was introduced in the 
context of thermodynamical systems having negative temperatures, where one can withdraw heat from a negative 
temperature reservoir and convert it completely to work, see e. g. [10], p. 176. Since thermodynamical systems 
with negatives masses of their components should have negative temperatures (e.g. [11]) there is a link between 
thermodynamical properties of such systems and the ones discussed in this paper. In particular, a runaway process 
occurring in a thermodynamical system having two subsystems containing particles with opposite signs of masses has 
been discussed in e.g. [12]. It has been mentioned that this process is analogous to the self-acceleration of a pair of 
particles with opposites signs of their masses. 
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